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Abstract. We present the analysis of a lOOksec INTEGRAL (3-100 keV) observation of the transient X-ray pulsar V 0332-H53 
in outburst. The source is pulsating at Ppuise = 4.3751 ± 0.0002 s with a clear double pulse from 6keV to 60keV. The average 
flux was ~550mCrab between 20keV and 60keV. We modeled the broad band continuum from 5keV to 100 keV with a 
power-law modified by an exponential cut off. We observe three cyclotron lines: the fundamental line at 24.9)° ] keV, the first 
ha rmonic at 50.5^°'! ke V as well as the second harmonic at 7L7)gg keV, thus confirming the discovery of the harmonic lines 
bv ICobum et alJ i200.5h in RXTE data. 
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1. Introduction 

V0332H-53 (RA=03‘'34“59)89, DEC=-h 53° 10’23”.6) is a re¬ 
current transient X-ra y pulsar. The system consists of the 
08-9 Ve star BQ Cam llNegueruela et al.lll999l:lBernacca et al.l 
1 1984 and a neutron star (NS) in an eccentri c (e - 0.31) 34.25 d 
orbit, spinning with a period of 4.375 sec jStella et alJll98-5ll . 
Since its discovery by the Vela 5B sat ellite during an outburst 
in 1973 llTerrell & Priedhorskylll984 . the system has so far 
exhibited three more outbursts in 1983 llTanakall 19831) . 1989 
(iMakinoi 19891) . and 2004. The peak luminosities during these 
outbursts were quite different: wh ile the 1973 outburst reached 
1.6 Crab (iNegueruela et al j 19991). the flux o f the 1983 outburst 
was about 10 times lower dUnger et al.|l992|) , and the 19 89 out¬ 
burst reached a flux of 0 .4 Crab )|Makishima_et_al .Il990l) . Based 
on a distance of 7 kpc jNegueruela et alJll999t) . which places 
the source beyond the Perseus spiral arm in the outer arm of 
the galaxy, the maximum X-ray luminosity reached during the 
1973 outburst was Lx > lO^^ergs^^ - close to the Eddington 
luminosity. 

Eo r the 1989 outburst observed by Ginga dMakishima et al.l 
1990ll . the X-ray spectrum has been modeled by a power-law 
plus a high energy cut off. The analysis further revealed a deep 
and broad absorption line at 28.5 keV which was interpreted as 
a cyclotron resonant scattering feature (CRSE). The 1989 data 
also revealed the presence of a quasi periodic oscillation (QPO) 


Send offprint requests to: 1. Kreykenbohm, 
e-mail: Ingo.Kreykenbohm @ obs.unige.ch 


with a centroid frequency of 0.05 Hz (iTakeshima et al.iri994 . 
The source furthermore e xhibits aperiodic s hort term variabil¬ 
ity similar to Cygnus X-1 dStella et al .119851) . 

The optical compani on BQ Cam was observed to b righten 
from January 2002 on bv iGoranskii & Barsukovald2004 . Since 
a similar brightening had been observed to precede the two pre¬ 
vious X-ray outbursts, they predicted an outburst of V 0332-1-53 
within the next one or two years. Indeed, the source entered 
an outburst in November 2004, after almost 15 years of qui¬ 
escen ce. recorded by the All-Sky Monitor (ASM. ISwank et all 
l2004 onboard the Rossi X-ray Timing Explorer (RXTE). The 
X-ray flux increased continuousl y and reached 1 & ab in the 2- 
12keV band in 2004 December dRemillardl2004 . Subsequent 
pointed observations with RXTE revealed the presence of 
three cyclot ron lines, at 26.34 + 0.03 keV, 49.1 + 0.2keV, and 
74 + 2keV dCobiirn et aLl200.5l) . 

In this Letter, we report on the 3-lOOkeV analysis based on 
the observation of V 0332-1-53 performed on 2005 January 6- 
10 by the INTEGRAL satellite. We describe the INTEGRAL 
observations in Sect. |2l Imaging, spectral, and timing analysis 
are presented in Sect.|3 and summarized in Sect.|4] 

2. INTEGRAL observations 

Although the peak of the outburst of V 0332-1-53 was reached 
in 2004 December, a dedicated target of opportunity (TOO) ob¬ 
servation could not be scheduled before 2005 January 5 due 
to satellite orientation constraints relative to the Sun. lOOksec 
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of good data were recorded for V 0332+53 between January 6 
and 10 - 30ksec in staring mode in revolution 272 (pointing 
75 on January 6th 08:47:45-18:14:00UT), 50ksec in hexago¬ 
nal mode in revolution 273 (pointings 66 to 82 from January 
8 at 22:18:12 to 9th at 15:45:47UT), and 20ksec in hexagonal 
mode in revolution 274 (pointings 4 to 9 on January 10 from 
03:20:30 to 09:22:11 UT). V 0332+53 was always in the field 
of view (FOV) of all three high energy instruments. 


3. Data analysis 

We analyzed data from all instruments on board INTEGRAL. 
The X-ray instrument JEM-X provides data in the 2-30 keV en¬ 
ergy band. The ISGRI layer of the imager IBIS covers the 20- 
800 keV band with an imaging resolution of 2 arcmin and about 
10% energy resolution at 20keV. The PICsIT layer provides 
data above 200 keV. The spectrometer SPI provides high en¬ 
ergy resolution (1 keV below lOOkeV) in the 20keV to 8 MeV 
energy range. The optical monitoring camera (OMC) provides 
the magnitude in the V-band for selected optical targets. 

Imaging, spectral, and timing analyses were performed us¬ 
ing version 4.2 of the Offline Scientific Analysis (OSA) pack¬ 
age as described in the cookbooks available with the OSA. 

3.1. Image analysis 

The analysis of V 0332+53 is facilitated by the fact that it is 
the only bright X-ray source in the FOV. The mean flux dur¬ 
ing the three observations is about 550mCrab in the 20-60 keV 
energy band in ISGRI and 700mCrab in the 3-lOkeV en¬ 
ergy band in JEM-X, consistent with RXTE/ASM. The flux 
at 20 keV during these observations is similar to th at measured 
in De cember 2004 during the Galactic Plane Scan dTiirler et all 
I 2 OO 4 I) . We note that, at lower energies, the RXTE/ASM has 
recorded a steady decline of the 3-lOkeV flux from a peak 
value of ~1 Crab at the end of 2004 December to ~0.6 Crab at 
the time of the INTEGRAL observation, the latter value being 
consistent with the JEM-X observation. 

The source was not detected by PICsIT; this is actually not 
surprising, knowing that the spectra of accreting X-ray pulsars 
are exponentially falling at higher ene rgies, making them diffl- 
cult to detect above 100keV (see, e.g.. lCoburn et al .l2002h . 

The OMC observed the optical companion of V 0332+53, 
BQ Cam at My ~ 15.4 + 0.2 consistent with simultaneous 
ground based observations dMasetti et al.ll200^ . The observa¬ 
tions confirm that BQ Cam is in an excited state. 


3.2. Spectral analysis 


We combined all data to obtain a spectrum with a very high 
signal to noise ratio. Spectra of accreting X-ray pulsars are, 
however, notoriously difficult to fit, especially if CRSFs are 
involved. Special care has also to be taken that no artificial 
residuals are introduced by the model (as is the case for, 
e.g., the broken power law and/or the high energy cut off, cf. 
Kreykenbohm et al. 19991) . We used the cutoffpl model of 


XSPEC which describes the overall shape of the continuum 
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Fig. 1. a Combined spectrum and model of data obtained with 
JEM-X (left), ISGRI (upper right), and SPI (lower right). Note 
that SPI is only shown for comparison, but is not used for the 
fits below. The dips of the cyclotron lines are already very ap¬ 
parent in the raw data. The continuum is modeled with the 
cutoffpl model, b model without any Gaussian lines applied. 
In c one Gaussian at 26.6/^] keV is included. In d another 
Gaussian is included at 47 .Uq 2 keV. e a third Gaussian is in¬ 
cluded 30.8 /;qj keV (see Table [0 resulting in a much better 
description of the spectrum between 20keV and 40keV, but 
significant deviations remain above 50keV. f the best fit using 
a fourth Gaussian at 71.7/;° ^ keV results in small residuals also 
at high energies. Note that all parameters mentioned previously 
are for the respective preliminary fits only. See Table^for the 
parameters of the final fit. 


reasonably well (see Fig.[0. We fit JEM-X and ISGRI data si¬ 
multaneously applying only a constant to allow for a free nor¬ 
malization between the two instruments. A systematic error of 
2% has been applied to account for the uncertainties of the re¬ 
sponse matrices of JEM-X and ISGRI. We did not use SPI as 
the fit is completely dominated by the ISGRI instrument above 
20keV anyway. For comparison, however, we also show the 
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Table 1. Fitted parameters from phase averaged spectra us¬ 
ing JEM-X and ISGRI spectra. All errors represent 90% confi¬ 
dence. 


Parameter 


2 CRSFs 

3 CRSFs 

r 


-0 46 

-0.20!“; 

Edit 

[keV] 

6 5 

- 0.1 

8.0!“ 

^Cyc.la 

[keV] 


94 Q +0-1 
^^•^-0.1 

Cr Cyc.la 

[keV] 

9 Q +0.1 
^•^-0.1 

9 4+0.1 

^•^-0.4 

rCyc.la 


1 n c -fO .01 
^ -0.01 

0 74 +0-^1 
'^•'^-0.02 

^Cyc,lb 

[keV] 

30.8!“ 

29.0!“ 

Cr Cyc.lb 

[keV] 

Q 7 -fO .2 
' -0.1 

'^•'^-0.1 

^Cyc,lb 


n 01 -f0.02 
-0.08 

1 9Q +0.01 
^•^°-o.oi 

^Cyc.2 

[keV] 

50 1 

JU.l _Q J 

50.5!“ 

c Cyc,2 

[keV] 

7 1 

' ‘ ^ -0.3 

8-112 

■rCyc,2 


1 CQ 4-0.03 
-0.01 

9 QQ +0.03 

_o .02 

Ecyc.S 

[keV] 

- 

71 1 ^^-^ 

' -0.8 

C Cyc,3 

[keV] 

- 

6 3 

- 0.4 

rcyc,3 


- 

1 Q 4-0.5 

1-0 n4 


(DOF) 

230.2 (147) 

170.4 (144) 


SPI data in Fig. [a, which shows the CRSFs at the same ener¬ 
gies as ISGRI. 

While the cutoffpl model hts the low energy data rea¬ 
sonably well without any photoelectric absorption, large ab¬ 
sorption line-like structures are present above 20keV. The ob¬ 
servation of harmonically spaced lines conhrms that these fea¬ 
tures originate from cyclotron resonant scattering. We use a 
multiplicative Gaussian absorption line at Ec =26.5;!;^ } keV 
to ht the strong fundamental CRSF betw een 20keV and 30keV 
(discovered by iMakishima et alJ[l990l). At ~50keV the h rst 
harmonic CRSF (discovered by RXTF. ICoburn et alJl2005l) is 
also extremely well detected (see Fig. QJ). After htting both 
CRSFs, many residuals still remain. This is due to t he non- 
Gauss ian shape of the funda mental, as observed by iMiharal 
( 1 19951) in the Ginga data and ICoburn et al.l ll2005l) in RXTE 
data. Using Lorentzian instead of Gaussian absorption lines 
does not improve the fit. To model the fundamental, we use 
a pair of Gaussians at 25.2 ;|;q[ keV and 30.8 ;^qj keV. The sig¬ 
nificant residuals visible at ~70keV (see Fig. m are identihed 
with the second harmonic CRSF (discovered by ICoburn et al.l 
I2OO5I) . They are fitted with a Gaussian at E - 71g keV (see 
Fig. m While the reported parameter values are for the re¬ 
spective preliminary hts, the hnal ht parameters are given in 
Table [2 

To further check the validity of our results, we used other 
spectral models to ht the data inclu ding the Negative Positive 
Exponential (NPEX. IMiharal 1 9^ . In all cases, we found that 
three CRSFs are present, although the actual line parameters 
depend on the shape of the continuum. We also note that the 
CRSFs observed at 25 keV and 50 keV are much more signif¬ 
icant than the deviations known to be present in the OSA 4.2 
response matrix or the energy calibration of ISGRI between 
20keV and 50keV; hence we are conhdent that the spectral be¬ 
havior is dominated by the source and not by calibration issues. 


100 
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Fig. 2. Power spectrum of the 50ksec observation in revolution 
273 of V 0332+53 (see Sect.EJ. The smaller peak at 0.22 Hz 
(4.375 s) is due to the pulse period while the higher peak at 
0.45 Hz (2.188 s) is due to the double-peaked pulse prohle of 
V 0332+53 (see Fig.|3. 

3.3. Timing anaiysis 

To determine the period of the pulsar we barycentered the event 
arrival times for ISGRI and JEM-X and derived power spectra 
with the xronos package for each of the three parts of the obser¬ 
vation (see Sect. El). No correction for orbital moti on has been 
applie d since the uncertainties of the ephemeris ilStella et all 
are too large to extrapolate over 20 years and the deriva¬ 
tion of a new ephemeris is out of the scope of this Letter. 

Two periodicities are clearly detected in the power spec¬ 
trum (see Fig. EJ: one peak at 0.2286 Hz associated with the 
spin of the NS and one at half period due to the double-peaked 
pulse prohle. The derived pulse periods are 4.3753 + 0.0001 s 
(Rev. 272), 4.37501 + 0.00006 s (Rev. 273), and 4.3749 + 
0.0002 s (Rev. 274). These variations are compatible with the 
effect expected from the orbital motion. 

We created pulse prohles in several energy bands. While the 
pulse prohle for Rev. 272 and 273 exhibits a clear double pulse 
at all energies, two behaviors are distinguishable in the pulse 
prohle for Rev. 274; below 6 keV only one pulse is detected, 
while two peaks are clearly observable at all energies above 
6 keV (see Fig. El- The secondary pulse is at a level of 60% 
to 75% amplitude of the main pulse (after subtracting the non- 
pulsed hux). At the same time, the pulsed fraction increases 
from 13% in 3-15 keV to 22% in 15-30keV and 35% in 30- 
60keV. 

We also performed rudimentary phase resolved spec¬ 
troscopy by deriving ISGRI spectra for the main pulse, the 
secondary pulse, and the pulse minimum using data from one 
pointing in Rev. 274. The normalized ratio between these spec¬ 
tra is consistent with unity between 20keV and 45keV, i.e., 
the spectrum in this energy range does not vary over the pulse. 
Above 45 keV the statistics are not sufficient from one point¬ 
ing. More detailed pulse phase resolved spectroscopy requires 
a valid ephemeris to correct for the binary motion. 
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Phase 


Fig. 3. Lightcurve of V 0332+53 in two JEM-X and two ISGRI 
energy bands in Rev. 274 (see Sect. |2j folded with Epuise = 
4.37492 s. Note the change of the pulse profile at low energies. 


4. Summary 


V 0332+53 is only the second X-ray source known so far that 
exhibits at least thr ee CRSFs; the othe r object, 4U 0115+63, 
exhibits hve lines llCoburn et al.ll20o3i . All three CRSFs in 
the INTEGRAL data are located close to those energies de¬ 
rived from RXTE. The fundamental CRSF has resolved struc¬ 
ture such that even two Gaussian lines are barely able to ht 
it. Comparing line shapes derived fr om pulse phase resolved 
spectra to Monte Carlo simulations of lArava & Hardinglill999l) 
would allow one to infer constraints on the geometry of the 
emission region. Such a study will be presented in a forthcom¬ 
ing paper. To derive the strength of the magnetic held, we need 
to derive the “true” centroid energy of the fundamental, which 
is difficult to obtain. The first harmonic at 50keV, however, 
is fit by a single Gaussian with small residuals (see Fig. 0. 
Assuming that the coupling factor between the fundamental 
and the first harmonic is indeed 2.0, we derive an energy for 
the fundamental of 25.3 + 0.1 keV, compatible with the cen¬ 
troid energy of the lower of the two Gaussians used to fit the 
25 keV feature. If we therefore assume that the energy of the 
fundamental CRSF is indeed 25.3 keV, that the NS has the 
canonical mass of 1.4 M©, and that the CRSFs originate close 
to surface of the NS we infer a strength of the magnetic held 
of B = (1 + z)/11.6 X Fcyc X lO'^G 2.7 x lO'^G with a 
gravitational redshift z of 0.25. The coupling factor between 
the fundamental CRSF and the second harmonic is 2.85 + 0.05, 
slightly smaller than the theoretical 3.0. The same has been 
observed for 4U 0115+63, where the shape of the fundamental 
line is also strongly distorted from that of a single Gaussian and 
the coupling factors are smaller than expected and is probably 


due to relativistic effects dHeindl et alJll99^ ISantangelo et alJ 
Il999ll . 

The derived pulse period is almost identical to the period 
of the 1983 outb urst (during which P - 4.3753 + 0.0001, 
Istella et al1ll98.5ll . indicating that the pulse period has, within 
the uncertainties arising from the lack of orbital correction, not 
changed significantly since 1983. 

The behavior of the source seems to be different below and 
above 6keV. The pulse profile exhibits two peaks at all en¬ 
ergies from 6 keV to 6 0keV. similar to other acc reting X-ray 
pulsars (e.g.. Vela X-1. iKrevkenbohm et alJl999h . while it ex¬ 
hibits only one peak below 6 keV in Rev. 274 and two peaks 
in Rev. 272 and 273. The behavior is different from the 1983 
outburst, where the sec ondary pulse was only present at very 
low energies (< 5 keV. I.Stella et alJll 98.51) . These authors also 
noticed a luminosity dependence of the pulse profile below 
lOkeV, two peaks being detected at higher luminosities and 
only one at lower luminosities. During the INTEGRAL ob¬ 
servation, the average luminosity in the 20-60 ke V band was 
~550m Crab - a high luminosity state according to lStella et al.l 
but we see a more complex behavior with one peak be¬ 
low 6 keV and two peaks above 6 keV in Rev. 274. The strength 
of the secondary peak is still 60% of the main peak between 
30keV and 60keV. Unfortunately, the statistics are not suf¬ 
ficient to determine a pulse profile at higher energy bands. 
Further analysis is in progress and will be reported in a forth¬ 
coming paper. 
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